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We review the recently developed critical minimum energy-subspace (CrMES) technique. This 
scheme produces an immense optimization of popular algorithms, such as the Wang-Landau (WL) 
and broad histogram methods, by predicting the essential part of the energy space necessary for 
the estimation of the critical behavior and provides a new route of critical exponent estimation. A 
powerful and efficient CrMES entropic sampling scheme is proposed as an alternative to the tra- 
ditional importance sampling methods. Utilizing the WL random walk process in the dominant 
energy subspace (CrMES- WL sampling) and using the WL approximation of the density of states 
and appropriate microcanonicai estimators we determine the magnetic properties of the 2D Ising 
model. Updating (E, M) histograms during the high level WL-iterations, we provide a comprehen- 
sive alternative scheme to the Metropolis algorithm and by applying this procedure we present a 
convincing analysis for the far tail regime of the order-parameter probability distribution. 
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I. INTRODUCTION 

Traditional Monte Carlo sampling methods have increased dramatically our understanding of the behavior of the 
standard classical statistical mechanics systems. The Metropolis method and its variants were, for many years, the 
main tools in condensed matter physics, particularly for the study of critical phenomena 0. [HQ. However, in many 
cases, such as for example in complex systems with effective complicated potentials, the Metropolis method and its 
variants are more or less inadequate methods 0, Q . Noteworthy, the traditional importance sampling methods are 
not very efficient in recording of the very small probabilities in the tails of the critical order-parameter distribution [j| . 
This is the reason why many different simulation approaches 0, 13 m recent times have failed in establishing the 
true tail behavior of this distribution. 

In Sec. [HI we briefly review a new simple technique, the CrMES method, which yields an immense speed up of 
all popular algorithms, which are used to sample the density of energy states (DOS) of a statistical system in the 
last decade. Provided that, the temperature range of interest is only the range around the critical temperature that 
determines all finite-size anomalies, our technique can be combined with any of the DOS methods [3. 13. M. fToL ITU Il2l 
Il3l HiL ITU Hil to speed up the simulations and establish a different methodology for the determination of the critical 
exponents by studying the finite-size scaling of the extensions of the dominant subspaces. In the same section we 
describe a CrMES- WL entropic sampling method and all magnetic properties are obtained by using the high-levels of 
the WL random walk process to determine appropriate microcanonicai estimators. The proposed method is combined 
with the N-fold way [l2l[la>[li| m order to improve efficiency and statistical reliability. Finally, in Sec. IIIII we study 
the order-parameter distribution of the 2D Ising model and we clarify its asymptotic behavior in the far tail regime. 
Our conclusions are summarized in Sec. IIVI 



II. THE CRMES WANG-LANDAU ENTROPIC SAMPLING SCHEME 

According to the CrMES method the total energy range (E min , E max ) is restricted using the definitions: 

{E-,E+), E±=E±A ± , A ± >0 (1) 
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with respect to the value E producing the maximum term in the partition function of the statistical model, for instance 
the Ising model, at some temperature of interest. Thus, the specific heat peaks are approximated by: 



C* L (E_,E+) = C2(A ± ) = N~ l T~' z { 



Z-^i^exp [$(£)]- [Z- 1 £ £exp [$(£?)] 



(2) 



= [S(E) - PE] - \S(E) -pE], Z = Y. ex P 



(3) 
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FIG. 1: Behavior of CrMMS for r = 10 -4 . The solid lines correspond to power law fits (a) for T = T c : a ~ 0.513(15) ■ L 1758(7) 
and (b) for T = T£ iX : S ~ 0.517(45) • L L741 ' 18> . The dotted line represents the law 3 = 0.5138 • L 1 ' 76 , which is very close to 
the fit 3~ 0.5145(90) • L 1 - 7497{50) corresponding to the average (H Tc + H T £ )/2. 

To implement the above restriction we request a specified accuracy by imposing the condition: 



C* T 



- 1 



< r 



(4) 



where r measures the relative error and will be set equal to a small number (r = 10~ 4 and/or r — 1CP 6 ), and C* L is 
the value of the maximum of the specific heat obtained by using the total energy range. With the help of a convenient 
definition pel l22j|. we can specify the minimum energy subspace satisfying the above condition. In fact, it has been 
numerically verified that the finite-size extensions (denoted by AE = min(£' + — E_)) close to a critical point obey 
the scaling law pfj| : 



(A£)^ 
L d 



(5) 



To locate the CrMES we may follow the method described in |2jj, or an even simpler restriction based on the 
energy probability density (Jti,c{E) oc &(E)). Using this later approach, we may define the end-points (E±) of the 
subspaces by simply comparing the corresponding probability densities with the maximum at the energy E: 



E±: exp{$(£; ± )} < r 



(6) 
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FIG. 2: Behavior of CrMMS for r = 1(T 6 . The solid lines correspond to power law fits (a) for T = T C : 3 ~ 0.528(17) ■ L 1761(7) 
and (b) for T = T£ x : 3 ~ 0.536(45) ■ L 1 - 745 ( 15 ). The dotted line represents the law 3 = 0.54 • L 1 ' 75 , which is very close to the 
fit 3 ~ 0.532(15) • L 1 - 7fi3(6) corresponding to the average (3 Tc + 3 T » )/2. 



Alternative definitions for the CrMES have been described and tested in Ref. [22| using the 2D Ising model 20, 22], 
the 3D Ising model HiHH and the Baxter- Wu model plli^. 

Let us now discuss the idea of producing accurate estimates for finite-size magnetic anomalies by using a simple 
method based on a WL random walk in an appropriately restricted energy subspace (E\,E2). Implementing this 
scheme we, at the same time, accumulate data for the two-parameter (E, M) histogram. A multi-range algorithm 0] 
is implemented to obtain the DOS and the (E, M) histograms in (E\, E%). The WL modification factor (fj) is reduced 

1/2 

at the jth iteration according to: f\ = e, fj — > fj_ 1 , j = 2 1 ...,Jfi n . The approximation of the DOS, in the last 
WL iteration, Gwl(E), and the high-level (j 3> 1) WL (E,M) histograms, Hwl(E, M), are used to estimate the 
magnetic properties in a temperature range, which is covered, by the restricted energy subspace (E\, E2) as: 

(Mn) E E (M n )EG(E)e-^ ^ E E e( El ,E 2) (M n )E,w L G WL (E)e^ 

The microcanonical averages (M n ) E are obtained from the Hwl(E, M) histograms as: 

(M n ) E S (M n )E,WL ^ r J? r? ' H WL (E) =Y j H wl {E,M) (8) 

and the summation in the magnetization M runs over all values generated in the restricted energy subspace (£1, E2). 

The accuracy of the magnetic properties obtained from the above averaging process depends on many factors. 
However, since the detailed balance condition depends on the control parameter {fj), we classify our recipes utilizing 
the j-range used for updating the (E, M) histogram during the WL process. The high-level recipes WL( Jj n ju Jfin) 
and their N-fold versions WL(N-fold: JN-fold, Jfin) give excellent estimates, as shown in detail in Ref. 22]. The 
extensions of the CrMES defined with the help of the thermal finite-size anomalies (specific heat maximum and 
energy cumulant minimum) , and also with the help of magnetic finite-size anomaly (susceptibility maximum) obey a 
clear logarithmic dependence on the lattice size [23, as should be expected for the 2D Ising model. 

Consider now the probability density of the order-parameter at some temperature of interest T: 



MM) = g -y> ™ ±2 

z2 E& ( El ,E 2 ) G WL(E)e f- 
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If M is the value that maximizes (0, we locate the critical minimum magnetic subspaces (CrMMS) by: 

~ P T * (M±) 

M± : L ' x _ < r (10) 
P T . x (M) 

and we should expect the following scaling law to hold: 

(AM)|, 

-J^LkLTi/" (11) 
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FIG. 3: Illustration of the universal scaling function p*(x) for L = 60 and L = 120. 

In Figs. lll2l we present the behavior of the extensions of the CrMMS, obtained using the magnetic space restriction 
(|1UH for two values of r, namely r = 10 -4 and r — 10~ 6 . These extensions were determined at the exact critical 
temperature T c and at the pseudocritical temperatures T£ of the susceptibility, as indicated in the figures. The 
power laws and the resulting exponents are given in the corresponding figure captions. It is clear that the law (|lll) is 
very well satisfied and the estimation of j/u via this route is quite good, giving ^jv = 1.7497(50) for r = 10~ 4 and 
7/V = 1.753(6) for r = 10~ 6 , respectively. It is of interest to point out that, an attempt to estimate j/v using this 
route and the Metropolis algorithm will most likely yield a marked underestimation, which will be more pronounced 
for the case r = 10~ 6 . This effect, has been discussed in Ref. [2^ and is a result of the very slow equilibration process 
of this algorithm in the far tail regime of the order-parameter distribution. Next, we consider the tail regime of 
the universal order-parameter distribution, we analyze its asymptotic behavior and finally examine the possibility of 
extracting estimates of the exponent S from this behavior. 




III. FAR TAIL REGIME OF THE ORDER-PARAMETER DISTRIBUTION 

Following Smith and Bruce we define the universal scaling form of the order-parameter density by: 

p* (x)dx ~ p(m)dm, x — vnj \J (m 2 ), m — M/N (12) 

Fig. shows p*{x) for L = 60 and L = 120 obtained via the WL(N-fold: 12-24) scheme. The tail regime of this 
universal distribution has been a matter of increasing interest in the last decade [H 0, 123L I24I |25| . The main theme is 
the verification of the following conjecture for the large- a; behavior of p*(x) |^,E|: 

p*{x) ~ Poo x^exp(-a 00 x s+1 ) (13) 



5 



with: 



*l> = S -^A ( 14 ) 



and pm, universal constants. The above hypothesis, its origin and its significance are discussed in some detail in 
Refs. EH. 
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FIG. 4: Behavior of estimates of the exponent ip for lattice sizes L = 60 — 140. The points correspond to the ratio of the best 
fit value of the exponent to the expected value 7 (to be compared with Fig. 7(b) of Ref. |g])- 



The studies of |5j, |6j have provided some evidence for this conjecture and in particular for the prefactor and the 
relation of the exponent ip to the critical exponent S. For the 2D Ising model the exponent should have the value 
ip — 7, if, of course, the prefactor hypothesis is valid. Smith and Bruce [5| provided numerical support for this value, 
but their study was not completely conclusive since it was carried out only for relatively small lattices (L — 32 and 
L = 64) and the x-window in which the value ip — 7 was observed was actually quite narrow. We now present 
results for several lattice sizes (L = 80, 100, 120, and L = 140) reinforcing this conjecture in a very wide x-window. 
Following Smith and Bruce we fix the exponent S in the exponential factor of Eq. I|13|l and fit our results (x > 1) 
in a;- windows, each one corresponding to 50 different magnetization values, sampled during the WL(N-fold: 12-24) 
process. Fig. 0] shows a very clear signature of the prefactor law Ijl4|l which upholds in a large x-window only for 
the large lattices (L = 120 and L = 140). On the other hand, for smaller lattice sizes (L = 60 — 100) the picture 
is similar to that presented in Ref. [5| and the expected value is obtained only in a small x-window. Finally, let us 
treat both exponents in the exponential and in the prefactor as free parameters, assuming however the validity of 
Eq. I|14|) . In this case we examine whether the far tail regime could be a possible route for an independent estimation 
of the exponent 6. Fig. [S] presents such an attempt for the case L = 120. The value obtained for ijj (and therefore 
for 6(= 2t/j + 1)) is accurate to the third decimal place [ip = 6.995(1.750)), despite the fact of the existing large error 
of the fit. In fact the shown x-window (x = 1.4 — 1.45) produces good estimates for all lattice sizes L — 60 — 140 as 
should be expected by thoroughly inspecting Fig. 0J 



IV. CONCLUDING REMARKS 



In the present work we described the basic ideas and formalism behind the recently developed by the au- 
thors [20L I2H |22| critical minimum energy(magnetization) subspace CrME(M)S technique, appearing as an alternative 
computational method for the estimation of critical behavior of statistical systems. We presume that this new promis- 
ing route will increase our comprehension of the development of the critical behavior as the lattice size increases and 
will further facilitate the estimation of critical exponents via finite-size scaling. It is hoped that the presented CrMES 
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FIG. 5: Independent estimation of the exponent S via fitting the universal scaling function p*{x) in a specific a;- window. 
The fitted line reads as p*(x) = p oa x i ' expl-a^x 2 ^ 2 ), yielding p x = 2.189(2.223), ip = 6.995(1.750) and aoo = 0.058(64), 
respectively. 

Wang-Landau entropic sampling will provide an efficient and reliable scheme for the study of complex systems, such 
as the random-field Ising model and systems with competing interactions, where the critical behavior is still contro- 
versial. Furthermore, the efficiency of our approach enabled us to present here reliable data for large lattices and to 
clarify the asymptotic tail behavior of the universal critical distribution of the order-parameter of the 2D Isingmodel. 
This should be contrasted to the fact that all previous studies of the far tail regime of this distribution yj |25j have 
considered rather small lattices. 
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